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We construct a general measure for the degree of non-Markovian behavior in open quantum 
systems. This measure is based on the trace distance which quantifies the distinguishability of 
quantum states. It represents a functional of the dynamical map describing the time evolution of 
physical states, and can be interpreted in terms of the information flow between the open system 
and its environment. The measure takes on nonzero values whenever there is a flow of information 
from the environment back to the open system, which is the key feature of non-Markovian dynamics. 
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The prototype of a Markov process in an open quan- 
tum system is given by a quantum dynamical semigroup, 
i.e., by the solutions of a master equation for the reduced 
density matrix with Lindblad structure [HE]- However, 
in realistic physical systems the assumption of a Marko- 
vian dynamics can only be an approximation that relies 
on a number of mostly rather drastic simplifications. In 
complex quantum systems one therefore often encoun- 
ters dynamical processes which deviate not only quanti- 
tatively but also qualitatively from the relatively simple 
behavior predicted by a Markovian time evolution [3J. 

In view of the large variety of conceptually different 
analytical methods and numerical simulation techniques 
that have been developed to treat non-Markovian sys- 
tems in recent years (see, for example, Refs. [31415]). the 
following questions arise: How can one rigorously define 
quantum non-Markovianity and how can one quantify 
the degree of non-Markovian behavior in a way which 
does not refer to any specific representation or approxi- 
mation of the dynamics, e.g. to a master equation with a 
given structure? In order to answer these questions one 
needs a measure for the non-Markovianity of the quan- 
tum dynamics of open systems which, in mathematical 
terms, represents a functional of the dynamical map that 
describes the time evolution of physical states. 

Here, we construct such a measure for non- 
Markovianity. This measure is based on the trace dis- 
tance of two quantum states which describes the prob- 
ability of successfully distinguishing these states. The 
basic idea underlying our construction is that Markovian 
processes tend to continuously reduce the distinguishabil- 
ity between any two states, while the essential property 
of non-Markovian behavior is the growth of this distin- 
guishability. Interpreting the loss of distinguishabilty of 
states as a flow of information from the open system to its 
environment, one is thus led to a simple, intuitive picture, 
namely that the key feature of non-Markovian dynamics 
is a reversed flow of information from the environment 
back to the open system. An important consequence of 



this picture is that the dynamical map of non-Markovian 
processes must necessarily be non-divisible, a property 
that is known to play also a decisive role in the classifi- 
cation of quantum channels [16 . 

To construct the measure for non-Markovianity we first 
need a measure for the distance of two quantum states 
pi and p2 ■ Such a measure is given by the trace distance 
(see, e.g., Ref. [IT]) which is defined by 

D{p 1 ,p 2 ) = ^tr|pi -p 2 \, (1) 

where \A\ = V A^A. The trace distance D represents 
a natural metric on the space of density matrices, i.e., 
on the space of physical states, satisfying < D < 1. 
Besides many other interesting properties, it has a clear 
physical interpretation in terms of the distinguishabil- 
ity of quantum states. Suppose that Alice prepares a 
quantum system in one of two states p\ and p 2 , each 
with probability | , and gives the system to Bob who per- 
forms a measurement to decide whether the system was 
in the state p\ or p 2 . One can show that the quantity 
\ [1 + .D(pi, P2)] is then equal to the probability that Bob 
can successfully identify the state of the system. Thus, 
the trace distance can be interpreted as a measure for the 
distinguishability of two quantum states. A further re- 
markable feature of the trace distance is given by the fact 
that all completely positive and trace preserving (CPT) 
maps $ are contractions for this metric j!8j . 

D($ Pl ,$p 2 )<D( Pl ,p 2 ). (2) 

This means that no trace preserving quantum operation 
can ever increase the distinguishability of two states. 

Suppose now that we have a quantum process given by 
a Markovian master equation, 

~At)=£p{t), (3) 
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with a generator in Lindblad form [TJ [5] , 



A iP A\ 



\{AAup} 



(4) 



involving a time-independent Hamiltonian H as well as 
time-independent Lindblad operators Ai and positive re- 
laxation rates ji > 0. Such a master equation leads to 
a dynamical semigroup of CPT maps $(<) = exp(Ct), 
t > 0, which describes the dynamics of the density ma- 
trix through the relation p(t) — $(t)p(0). By use of the 
semigroup property $(t + £) = $(r)$(t) it easily follows 
from Eq. ^ that for all r, t > we have 



D( Pi (t + t),p 2 {r + t)) < D{ Pl (t),p 2 (t)), 



(5) 



where pipit) = <&(t) P i^{0). Thus, for all quantum dy- 
namical semigroups the trace distance of the states 
Pifl{t), corresponding to any fixed pair of initial states 
Pi,2(0), is a monotonically decreasing function of time. 
This is a general feature of quantum Markov processes, 
implying that under a Markovian evolution any two ini- 
tial states generally become less and less distinguishable 
as time increases. We can interpret this loss of distin- 
guishability as a certain flow of information from the 
system to the environment which continuously reduces 
our ability to distinguish the given states. 

The inequality ^ holds for a much larger class of 
quantum processes than those described by a master 
equation of the form ([3]). In fact, suppose we have a 
time-local master equation of the form 



dt 



p(t) = K{t)p{t) 



(6) 



with a time-dependent generator K.(t). One can show 
that in order to preserve the Hermiticity and trace of the 
density matrix this generator must be of the form [TJ [7] 



K{t)p = -i[H(t),p] 



(7) 



$>(*) \Mt) P Al(t)-^{Al(t)Mt),p} 



where the Hamiltonian H(t), the Lindblad operators 
Ai[t) and the relaxation rates "/i(t) depend on time. If 
the relaxation rates are positive functions, "fi(t) > 0, the 
generator is seen to be in the Lindblad form Q for 
each fixed t > 0. Such processes with ji(t) > may 
be called time-dependent Markovian although the corre- 
sponding dynamical maps do not lead to a quantum 
dynamical semigroup. With the help of the chronological 
time-ordering operator T we can define a two-parameter 
family of CPT maps $(^2,^1) by means of 



*(i 2 ,ti) =Texp 



dt'fC(t') 



ti 



(8) 



The dynamical map which transforms the initial states 
at time into the states at time t can then be written as 



= 4>(t , 0) . The important point to note is that this 
dynamical map has the property of being divisible in the 
sense that for all r, t > the CPT map $(r + t, 0) can be 
written as composition of the two CPT maps $(r + t, t) 
and $(t,0), 



$(T + t,0) = *(T + f,t)$(t,0). 



(9) 



We remark that for a dynamical semigroup on has 
$(*2j*i) = $(h - h) sucn that Eq. ^ reduces to 
$(t+*) = $(r)$(t). Since in Eq. ^ not only $(t+£,G) 
and <f>(t, 0) but also $(r+t, t) is a CPT map, we conclude 
that the relation ([5| holds true for all time-dependent 
Markovian quantum processes defined by the master 
equation ([6| with ji(t) > 0. 

We define the rate of change of the trace distance by 



cr(t, Pl , 2 (0)) = j t D( Pl (t),p 2 (t)) 



(10) 



For a given quantum process $(t), this quantity depends 
on time t and on the initial states pi^O) with corre- 
sponding time-evolutions pxzit) — ®(t, 0)pi,2(0). As has 
been demonstrated above, we have a < for all quantum 
processes for which the divisibility property ^ holds, 
i.e., for all dynamical semigroups and all time-dependent 
Markovian processes. We remark that Eq. ^ not only 
holds for CPT maps, but also for the larger class of pos- 
itive and trace-preserving maps |18| . Thus, a < holds 
true also for Markovian master equations which are not 
in Lindlbad form but preserve positivity. 

There are however many physical processes for which 
a is larger than zero for certain times. It is this type of 
processes which we define as non-Markovian. Hence, a 
process is said to be non-Markovian if there exists a pair 
of initial states /91.2(C)) and a certain time t such that 
Pi, 2(0)) > 0. Physically, this means that for non- 
Markovian dynamics the distinguishability of the pair 
of states increases at certain times. We interpret this 
as a flow of information from the environment back to 
the system which enhances the possibility of distinguish- 
ing the two states. While Markovian processes tend to 
wash out more and more characteristic features of the 
two states, non-Markovian processes lead to an uncover- 
ing of these features. We emphasize that the temporary 
backflow of information represents a natural feature oc- 
curring in many physical systems which does not imply 
that there is no thermalization for long times. 

How can one construct a measure for non-Markovianity 
on the basis of this definition? Clearly, such a quantity 
should measure the total increase of the distinguishability 
over the whole time-evolution, i.e., the total amount of 
information which flows from the environment back to 
the system. This suggests defining a measure J\f(&) for 
the non-Markovianity of the quantum process $(t) by 
means of the relation 



J\f(&) = max 



dta(t,p 1<2 (0)). 



(11) 



er>0 
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Here, the time-integration is extended over all time in- 
tervals (a,,6j) in which a is positive, and the maximum 



is taken over all pairs of initial states. In view of Eq. ( 10 ) 
we can thus write this definition as 

= max V \D(pi(bi),p2(bi))-D(pi(ai),p 2 (a,i)) 
01,2(0) L 

(12) 

To calculate this quantity one first determines for any 
pair of initial states the total growth of the trace dis- 
tance over each time interval (a,i,bi) and sums up the 
contributions of all intervals. Af(&) is then obtained by 
determining the maximum over all pairs of initial states. 





FIG. 1: (Color online) (a) The rate of change a of the trace 
distance as a function of time t and detuning A for the initial 
pair of states p x (0) = |+)(+| and p 2 (0) = |-)(-|. (b) The 
black regions represent the regions in which a is positive. 



By construction, we have Af($) = for all processes 
which have the divisibility property (|9| . In the following 
we discuss two simple examples for which our measure 
of non-Markovianity is greater than zero. The aim is to 
illustrate how to determine this quantity in specific cases 
and how non-Markovianity is related to a violation of 
the divisibility property and to the emergence of negative 
rates in master equations of the structure ^ . 

The first example describes a two-level system with 
excited state |+) and ground state |— ) which interacts 
with a reservoir of field modes. The exact interaction 
picture master equation [3] describing the dynamics of 
the density matrix is of form of Eq. ^ with the gen- 
erator Q, where H(t) = and we have only a sin- 
gle time-independent Lindblad operator A = <r_ and 
a time-dependent rate 7(t). The function j(t) is de- 
termined by the spectral density J(oS) of the reservoir. 
We investigate the case of a Lorentzian spectral density 
J(uS) = 7oA 2 /27r[(w — A - uj) 2 + A 2 ], the center of which 
is detuned from the transition frequency loq of the two- 
level system by an amount A, and work in the weak 
coupling limit 7o/A = 0.01 (damped Jaynes-Cummings 
model). For sufficiently large detunings A, the function 
7(f) then describes an exponentially damped oscillation 
and takes on negative values within certain intervals of 
time corresponding to a revival of the coherence in the 
system [3J 2] . We emphasize that this does not imply a 
violation of the complete positivity of the corresponding 
dynamical map because the necessary and sufficient 



condition for the complete positivity of Q(t) is given by 
r(i) = J Q dt'f(t') > 0, which is indeed satisfied here. 

However, the trace distance increases for those t for 
which 7(t) < 0, i.e., we have a(t, pi,2(0)) > for these 
times. This point is illustrated in Fig. [I] which shows a as 
a function of time t and detuning A for the pure initial 
states pi(0) = l+K+l and p 2 {0) = |-)(-|. For these 
initial states one finds the simple expression 



ff(t,Pi,2(Q)) = -7(t)exp[-r(t)], 



(13) 



which shows that a positive a and an increase of the 
trace distance is linked to a negative rate in the master 
equation. Thus, the appearance of negative rates signifies 
a violation of the divisibility property ^ and a flow of 
information from the environment back to the system. 
The maximization over the pair of initial states pi.2(0) 



in expression (111 can be performed by drawing a suf- 



ficiently large sample of random pairs of initial states. 
The results are shown in Fig. [2] and provide strong nu- 
merical evidence that the maximum is attained for the 
initial states pi(0) = |+)(+| and p 2 (0) = \-)(-\. This 
result could have been expected since p2(0) represents 
the invariant state and pi(0) has the largest distance to 
this state. According to Fig. [2] J\f($) exhibits a non- 
monotonic behavior: The increase of the number of in- 
tervals in which a > is overcompensated for large A by 
the decrease of the size of a in these intervals. 
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FIG. 2: (Color online) The non-Markovianity _/V($) for the 
damped Janes-Cummings model as a function of the detuning 
A. Blue dots: 1000 randomly drawn pairs of pure and mixed 
initial states. Red circles: The initial pair pi(0) = |+)(+| and 
/°2(0) = |— )(— | which leads to the maximum in Eq. \Ll\ . 



For the previous example the non-Markovianity J\f(&) 
was found to take on finite, positive values. Our second 
example represents a rather extreme case, demonstrating 
that there are also processes for which Af(&) is infinite. 
We consider a central spin with Pauli operator <x which 
interacts with a bath of N spins with Pauli operators er( fc ) 
through the Hamiltonian H = A^2^ =1 u z ai k \ where A 
is a coupling constant. This simple model can easily be 
solved exactly [T!5] ■ Assuming the initial state of the bath 
to be a maximally mixed state, one finds that the pop- 
ulations of the density matrix of the central spin stay 
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constant in time, while the coherences are multiplied by 
the factor /(f) = cos^^Af). This leads to a simple for- 
mula for the trace distance of the states pi(t) and P 2(t), 



D( Pl (t),p 2 (t)) = ^/a2 + P(t)\b\< 



(14) 



where a = pf + (0) — p% (0) denotes the difference of the 
populations, and b = p\ (0) — p\~ (0) the difference of 
the coherences of the two initial states. It follows that 
the trace distance oscillates periodically between -D ma x = 
y/a 2 + \b\ 2 and -D m in = \a\. This can be interpreted as 
a periodic oscillation of the distinguishability of the two 
states, i.e., as a periodic exchange of information between 
the central spin and the spin bath. 

The maximal growth of the trace distance occurs if one 
takes as initial states the two eigenstates of a x , or any 
other pair of states corresponding to antipodal points on 
the equator of the Bloch sphere, such that a = and 
\b\ = 1. The trace distance then oscillates periodically 



between the values 1 and 0. The sum in Eq. ( 12 ) therefore 
diverges and we obtain A/"($) = +oo, which implies that 
a Markovian approximation of the system dynamics is 
never possible. One can formally write a master equation 
of the form Q with H = 0, a single Lindblad operator 
A = o z and the rate 7(f) = AN tan(2 At), which shows 
again the connection between the growth of the trace 
distance and the emergence of negative rates. 

Summarizing, we have constructed a measure A/"($) for 
the non-Markovianity of quantum processes in open sys- 
tems. The definition (111 of the measure neither relies on 



any specific representation or approximation of the dy- 
namics, nor does it presuppose the existence of a master 
equation or of invariant states. The exact determina- 
tion of the measure generally requires solving the com- 
plete reduced system dynamics which could be a difficult 
task for more complex problems. However, any observed 
growth of the trace distance is a clear signature for non- 
Markovian behavior and leads to a lower bound for J\f(&). 
The measure for non-Markovianity introduced here could 
therefore be useful also for the experimental validation 
of theoretical models or approximation schemes. To de- 
tect non-Markovianity experimentally one has to perform 
a state tomography on different ensembles at different 
times in order to decide whether or not the trace dis- 
tance has increased. A great advantage of the present 
approach is given by the fact that it allows to plan exper- 
iments which test non-Markovianianity without know- 



ing anything about the properties of the environment or 
about the structure of the system-environment interac- 
tion. Hence, we think that our measure is a useful tool 
for the characterization of non-Markovianity, both in the- 
oretical descriptions and in experiments. 
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